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^ ■ Abstract 

All the possible schemes of neutrino mixing with four massive neutrinos 
^ ■ inspired by the existing experimental indications in favor of neutrino 

mixing are considered in a model independent way. Assuming that in 
short-baseline experiments only one mass-squared difference is relevant, 
it is shown that the scheme with a neutrino mass hierarchy is not com- 
patible with the experimental results. Only two schemes with two pairs 
of neutrinos with close masses separated by a mass difference of the 
order of 1 eV are in agreement with the results of all experiments. One 
of these schemes leads to possibly observable effects in and {f3(3)oi, 
experiments. 



1 Introduction 



Neutrino masses and neutrino mixing are natural and plausible phenomena of modern 
gauge theories (see, for example, Ref.|^). However, for the time being, the values of 
the neutrino masses and mixing angles cannot be predicted on theoretical grounds. The 
determination of these quantities is the key problem of today's experimental neutrino 
physics. 

At present there are several indications in favor of neutrino masses and mixing. One 
of the most important indications comes from solar neutrino experiments (Homestake 0, 
Kamiokande 0, GALLEX and SAGE ||^). As it is well-known, in all four presently 
operating solar neutrino experiments the observed event rates are significantly smaller 
than the values predicted by the Standard Solar Model (SSM) 0. Moreover, if the sur- 
vival probability of solar u^s is equal to one, the data of different experiments cannot be 
explained even if the total neutrino fluxes are considered as free parameters 0. Assuming 
the validity of the SSM, the experimental data can be described by the MSW matter ef- 
fect IP] for Am^ ~ 10^^ eV or by vacuum oscillations in the case of Am^ ~ 10^^'' eV 
||T0| (Am^ is the neutrino mass-squared difference). 

The second indication in favor of neutrino mixing comes from the data of the Ka- 
miokande 1MB |]l2l and Soudan JTsf atmospheric neutrino experiments. The ratio 



of muon-like to electron-like events measured in these experiments is less than the ex- 
pected ratio. The sub-GeV and multi-GeV data of the Kamiokande collaboration can be 
explained by ^ u^- or ^ oscillations with Am^ ~ 10~^eV^. 

Finally, in the LSND experiment |]14| u^p —>■ n events produced by neutrinos orig- 



inating from vr"*" and subsequent /i"*" decays at rest were observed. These events can be 
explained by ^ v^, oscillations with Am^ ~ 1 eV^ . 

Therefore, from the existing experimental data we have three different indications in 
favor of non-zero neutrino masses, each with a particular scale of Am^. However, we 
must also take into account the fact that in several short-baseline experiments neutrino 
oscillations were not observed. The resuls of these experiments allow to exclude large 
regions in the space of the neutrino oscillation parameters. 

In the present paper we will examine what information on the neutrino mass spectrum 
can be inferred from the results of all short-baseline neutrino oscillation experiments if we 
also take into account the results of the solar and atmospheric neutrino experiments. We 
will present a general discussion which does not assume a definite model of neutrino mix- 
ing. We will show that the experimental results favor two rather particular possibilities 
for the neutrino mass spectrum. 

We will start with the presentation of the general formulas for the neutrino transition 
probabilities in short-baseline neutrino oscillation experiments (for details see Ref . [p!5[] ) . 
Our basic assumption is that only one neutrino mass-squared difference is relevant for 
short-baseline neutrino oscillations. This assumption means that the neutrino mass spec- 
trum consists of two groups of close masses, separated by a mass difference in the eV 
range. Denoting the neutrinos of the two groups by i^i, . . . , z/,, and z/r+i, ... , z^n, respec- 
tively, with masses mi < . . . and m^+i < . . . < m„, we can quantify our assumption 
by 

Am,^iL ^ , Aml:L ^ , , 

— < 1 for %<r and ^ < 1 for i > r + 1 , (1) 

2p 2p 
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where Am?- = mj — m?, L is the distance between the neutrino source and detector and 

I— K 

p is the neutrino momentum. We would hke to emphasize that the inequahties (|l]) are 
well satisfied for short-baseline experiments with the mass squared differences relevant 
for the explanation of the solar and atmospheric neutrino experiments. 

Under the assumption of validity of the inequalities (|1]) , the amplitude of the transition 
is given by 



A, 



i>r+l 



exp —t 



Am^L 



2p 



(2) 



Here Am^ = — ml, U is the unitary n x n mixing matrix and u^, are any active 
or sterile neutrinos. (Note that the number of sterile states is n — 3.) From Eq.(0), for 
the probability of the transition z/^ — z//^ (a 7^ [3) we obtain 



1 , / Am^L 

Pv^^v^ = -Aa-/3 I 1 - COS 



2 

where the oscillation amplitude Aa-f3 is given by 



2p 



(3) 
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i>r+l 



A 



(4) 



The survival probability of is calculated as 



1 - ^Bn 



1 — COS ■ 



Am^L 



2p 



(5) 



where 



B. 



i>r+l \ i>r+l ) 



(6) 



Using the unitarity of the mixing matrix, from Eqs.(HD and (BI) we have also 



A 



5. 



4Eit^mr |i-Ei^-n • 

i<r \ i<r / 



(7) 



The expressions describe the transitions between all possible neutrino states, 

whether active or sterile. Let us stress that in the scheme under consideration the oscil- 
lations in all channels are characterized by the same oscillation length Lose = 47rp/Am^. 
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2 Schemes with three massive neutrinos 



Neutrino oscillations in the scheme with three massive neutrinos and a mass hierarchy 
were considered in many papers |T^, ^7\, |I8|. In Ref.|r7| it was shown that the results 
of the LSND experiment are compatible with the data of all the other experiments on 
the search for neutrino oscillations and the data of solar neutrino experiments only if the 
element I^T^sl of the neutrino mixing matrix is large (close to one) and jf/esl; l^rsl are 
small. 

In Refs.||18|| the possibility was considered that the same Am^ ~ 0.3 eV^ is relevant for 
the LSND effect and for the atmospheric neutrino anomaly. This possibility seems to be 
rather marginal. In fact, from the results of the LSND and Bugey [|I^ experiments 
it follows that the transition amplitude A^-e is less than 4 x 10~^. Therefore, in order to 
explain the sub-GeV Kamiokande data one needs a large amplitude A^.^ ~ -B^;^. Taking 
into account the limits for 5^.^ that follow from the CDHS disappearance experiment 
one finds Am^ < 0.4 eV^. On the other hand, the results of the LSND and Bugey 



experiments require Am > 0.3 eV . Thus, the results of all short-baseline oscillation 
experiments are compatible with the Kamiokande sub-GeV data only if the value of Am^ 
lies in the very narrow interval 0.3 eV^ < Am^ < 0.4 eV^ and A^.,, ~ 3 x 10~^. At such 
large values of Am^ the cosine in the expression (H) practically disappears for atmospheric 
neutrinos because of the averaging over energy and distance. This would mean that there 
cannot be an zenithal dependence of the double ratio R = (;u/e)data/(/^/e)MC of atmo- 
spheric muon and electron events ((/i/e)Mc is the Monte-Carlo calculated ratio of muon 
and electron events without neutrino oscillations). However, some zenithal dependence 
of this double ratio is indicated by the multi-GeV Kamiokande data. 

The existence of an zenithal dependence of the double ratio R will be checked soon 
by the on-going Super-Kamiokande experiment ||2l|. Let us also mention that ^ Ur 



oscillations with large amplitude and Am ~ 0.3 eV could be tested by CHORUS p2 



NOMAD [H and especially by the COSMOS H experiment. 

Up to now we have discussed a scheme with three neutrinos and a mass hierarchy. 
The other possible scheme with three neutrinos is a scheme where the spectrum has the 
form mi <^ 1712 < ||2^. With the same arguments as above it is possible to show |T5| 
that this scheme is also disfavored by the experimental results. 



3 Four massive neutrinos with a mass hierarchy 

We will now assume that the anomaly in the sub-GeV and multi-GeV Kamiokande atmo- 
spheric neutrino data is due to neutrino oscillations, for which a scale of Am^ ~ 10~^ eV^ 



is required |]TT|. In order to accomodate also the solar neutrino data and the LSND data, 
it is necessary to assume that the flavor neutrino fields are superpositions of (at least) four 
massive neutrino fields. We will consider now all possible schemes with four neutrinos 
and the dominance of one Am^ in short-baseline experiments. We will start with the case 
of a neutrino mass hierarchy, mi <^ m2 <^ ^ m^, with Am2i and Am|2 relevant for 
the suppression of the flux of solar neutrinos and for the atmospheric neutrino anomaly, 
respectively. This case corresponds to = 4 and r = 3 in the formulas (|3)-(0) and the 
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oscillation amplitudes are given by 



Ao,.,p = A\Up^\^\U^^W (9) 

Sa;„=4|f/«4|'(l-|t/a4|') . (10) 

We will consider the ranee 0.3 eV^ < Am^ < lO^eV^ which covers the sensitivity of 
all short-baseline experiments. At any fixed value of Am^, from the exclusion plots of 
the Bugey |jl9|, CDHS and CCFR |^ disappearance experiments we have 



Sa;a<<„ (« = 6, /i) . (11) 

The values of B^.^ and 5°.^ can be obtained from the corresponding exclusion curves. 
From Eqs.([10|) and (|lT]) we find that the elements |^7o4p must satisfy one of the two 
inequalities 

\Uo,A? <al or |?7«4r > 1 - a° (a = e, /x) , (12) 
where (see Ref. ]r7|) 

«° = ^ (i - \l^-Ko) ■ (13) 

In the range of Am^ considered here Og and a° are small (Og < 4 x 10^^, a° < 10^^). We 
will show now that solar neutrino data and the atmospheric neutrino anomaly exclude 
large values of |f/e4p and |f/^4p. In fact, the average probability of solar neutrinos to 



survive is given by (see Refs.[p7 



= E + f 1 - E \uA ^iSe , (14) 

j=3,4 V j=3,4 / 



where Pu^lJ^^ is the survival probability due to the mixing of with vi and V2- If 
|f^e4p > 1 ~ ^e, from ( p!4D we have P®^^^ > 0.92 for all solar neutrino energies. Such a 
large lower bound is not compatible with the solar neutrino data. 
The average probability of atmospheric i/^'s to survive is given by 

where Pul'^tl is the z/^ survival probability^ due to the mixing of z/^ with z/3 and 1^2, z^i- 
The double ratio of atmospheric muon and electron events is given by 

patm I patm 
Datm _|_ J- Patm ' ^ ' 



1 n(l,2;3) • • , 

P>)^~^ul IS given by 



'^u^'^u /ITT I9.ITT I9.ITT i9\2l 



U,i\' + \U,2\' + \U,,\^Y V 2p 
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where r is the ratio of muon and electron events calculated without neutrino oscillations. 
For the Kamiokande sub-GeV events r ~ 1.57 and there is no zenithal dependence of the 
double ratio R |Tl|. This means that the oscillatory terms in the probabilities P^*™^^ 
= e,/i) in Eg . ([T6|) disappear because of the averaging over energy and distance. In 
this case P,?*™„ = -P,?*™„ < 1 — P,?*™„ . From Eq.(O), we obtain for R the lower bound 

Let us consider now the case |f/^4p > 1 — a°, which implies 

^C^.. = ■ (18) 

The solid curve in Fig.|I| depicts the corresponding value of i?™™ for the Kamiokande sub- 
GeV events as a function of Am^. The shadowed horizontal band in Fig.|I] represents the 
90% CL limits of the double ratio R for the sub-GeV events determined in the Kamiokande 
experiment [^. From Fig.|I] one can see that, for Am^ > 0.4 eV^ and |f/^4p > 1 — 
the value of i?™"^ is bigger than the experimental upper limit for R. 

With the argument presented above we cannot exclude large values of |f//i4p in the 
small interval 0.3 eV^ < Am^ < 0.4 eV^. However, if the atmospheric neutrino anomaly 
is mainly due to a deficit of /i-like events (as indicated by a comparison of the experi- 



mental data |11, 12, 13] with the existing calculations of the atmosperic neutrino fluxes 



p8| , with the exception of the calculation presented in Ref. p9|), large values of |t/^4p are 



incompatible with the data for all values of Am . This is connected with the fact that 
for large values of |?7^4p the coefficient (1 — |?7^4p)^ of the survival probability Pl^'^ll 
in Eq.(|l^) is very small and it is not possible to explain the zenithal dependence ob- 
served in the Kamiokande multi-GeV data ||ll|] (for example, at Am^ = 0.3 eV^ we have 
(l-|^M4p)'<0.06). 

Thus, using the results of the reactor and accelerator disappearance experiments and 
taking into account the solar neutrino data and the atmospheric neutrino anomaly, we 
come to the conclusion that 

\Ue4\^ < a° and \U^i\^ < . (19) 

We will consider now z/^ ^ oscillations. From Eqs.(^ and (W) we have 



A^-e = 4\Ue4rK4r < . (20) 

Thus, the upper bound for the amphtude A^.g is quadratic in the small quantities a^, 



aj;, and ^ oscillations must be strongly suppressed if there is a neutrino mass 
hierarchy. In Fig.^ the limit ( ^0|) is presented as the curve passing through the circles. 
The 90% CL exclusion regions found in the Ue disappearance Bugey experiment and in 
the Ufj^ — > Ue appearance BNL E776 |3^ and KARMEN |^ experiments are limited in 



Fig.0 by the dashed, dot-dashed and dot-dot-dashed curves, respectively. The shadowed 
region in Fig.|^ is the region of the parameters Am^ and A^.g which is allowed at 90% 
CL by the LSND experiment. It is seen from Fig.^ that the region allowed by LSND is 
inside of the regions that are forbidden by the results of all the other experiments. Thus, 
we come to the conclusion that a mass hierarchy of four neutrinos is not compatible with 
the results of all neutrino oscillation experiments. 
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4 Four massive neutrinos with non-hierarchial mass 
spectra 

If the neutrino masses satisfy the inequahties 

mi<^m2<m3<mi, (21) 

with Am32 and Am^g relevant for the suppression of the solar u^s and for the atmospheric 
neutrino anomaly, respectively, the short-baseline oscillation amplitudes are given by 
Eqs.(|^) and (|TD|) with the change lUail"^ — ^ Arguments similar to those presented 

in Section |^ lead us to the conclusion that the mass spectrum ( pT]) is disfavored by the 
experimental data. In a similar manner one can demonstrate that all possible schemes 
with mass spectra in which three masses are clustered and one mass is separated from 
the cluster by the ~ 1 eV gap needed for the explanation of the LSND data are not 
compatible with the results of all neutrino oscillation experiments. 

Now we are left only with two possible neutrino mass spectra in which the four 
neutrino masses appear in two pairs separated by ~ 1 eV: 

atm solar solar atm 

(A) nil < ^2 ^ m3 < ml and (B) mi < m-l ^ m^ < ml . (22) 
, ^ ' ^ ^ ' 

LSND LSND 

The possible effects of these neutrino mass spectra have been discussed in Refs.[^ |T5[. 
We will show now that the schemes with the spectra (A) and (B) are compatible with 
the results of all neutrino oscillation experiments. Let us define the quantities 

c„ = ^ (a = e,/i). (23) 

i=l,2 

For both schemes (A) and (B) the amplitude Ba-^a is given by (see Eq.(^ with n = 4 and 
r = 2) 

B^,^ = 4c„(l - Ca) . (24) 

From the results of reactor and accelerator disappearance experiments it follows that the 
parameters Cq, must satisfy one of the two inequalities 

Ca < a° or Ca>l- a° (a = e, (x) , (25) 



where is given by Eg. ([131) . 

Let us first consider the scheme (A). For the probabilities of solar u^s and atmospheric 



i/^'s to survive we have 



^.%.= Elf^-l' + (l-^^)'^i'^l (26) 



1=1,2 

and 



P:%,. = (1 - + cJP£i!.l^ . (27) 
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If Ce > 1 — ctg, from Eq.(p6D it follows that the survival probability of solar u^s, 
practically does not depend on the neutrino energy and 



P®^^ >0.5 



(28) 



This is disfavored by the solar neutrino data |^3|. If < aj^, from Eq. 

, 2 



we have 



(29) 



which is not compatible with the Kamiokande atmospheric neutrino data for the reasons 
discussed in Section |^. Thus, in order to accomodate the solar neutrino data and the 
atmospheric neutrino anomaly, from the four possibilities (^) we must choose 



c. < 



and 



(30) 



Let us consider now z/„ ^ z/g oscillations. Using the Cauchy-Schwarz inequality, from 



Eq.(H|) (with n = A and r = 2) and Eq.(E3|), for both schemes (A) and (B) we find 



4 = 1,2 



(31) 



From Eqs.(|OD and (|3T|) it follows that the upper bound for A^-e is linear in the small 
quantity a^. Since ctg ^ 5 x 10~^ for all values of Am^, in the case of the scheme (A) the 
limit (0) is compatible with the results of the LSND experiment. 

In the case of the scheme (B), the solar neutrino problem and the atmospheric neutrino 
anomaly can be explained by neutrino oscillations only if 



Cp > 1 - 



and 



(32) 



From Eqs.(pl|) and (^) it follows that the scheme (B) is also compatible with the results 
of the LSND experiment. 

The schemes (A) and (B) lead to different consequences for the experiments on the 
mesurement of the neutrino mass through the investigation of the end-point part of the 

/5-spectrum and for the experiments on the search for neutrinoless double-/? decay 
((/3/3)o;y). In fact, for the whole range of Am^ considered here we have 



(A) 5^|t/e.P>l-a° 

i=3,4 

(B) Yl ^ «° ■ 



(33) 
(34) 



i=3,4 



From Eq.(|33|) it follows that in the case of the scheme (A) the neutrino mass that enters 
in the usual expression for the j3 spectrum of decay (see Ref.||3^) is approximately 
equal to the "LSND mass" 1714^: 



my('^H) ~ m4 . 



(35) 



If the scheme (B) is realized in nature and mi, m2 are very small, the mass measured 
in experiments is at least two order of magnitude smaller than m^. 
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If massive neutrinos are Majorana particles, {PP)ou decay is possible. In the scheme 



(A), the effective neutrino "mass" |(m) 
is given by 

l{n>>l 



1=1 



i=3,4 



m4 



that is measured in {PP)ou decay 



(36) 



We have 



|(m)| ~ m/^sjl 



4|f/e4P(l 



sm 



(37) 



where is the difference of the phases of f/es and UeA- Depending on the value of the 
phase 0, the quantity | (m) | has a value in the range 



\2\U, 



e4| 



l| nii < \ {m) \ < . 



(3^ 



The upper and lower bounds in Eq.(PBD correspond, respectively, to the cases of equal 
and opposite CP parities of z/3 and z/4 (see Ref.[|15|). From Eq.(|34D it follows that, if mi 
and m2 are very small, in the case of the mass spectrum (B) the expected value of | (m) | 
is very small. 

Thus, the experiments on the investigation of the effects of neutrino masses with the 
measurement of the end-point part of the /3-spectrum of and with the search for {l3(3)Qy 
decay could allow to distinguish between the possibilities (A) or (B) for the neutrino mass 
spectrum. 



5 Conclusions 

In this paper we have discussed the possible form of the neutrino mass spectrum that can 
be inferred from the results of all neutrino oscillation experiments, including the solar 
and atmospheric neutrino experiments. In this investigation we assumed only that one 
neutrino mass squared difference is relevant in short-baseline oscillation experiments. We 
have argued that it is unlikely that all data can be fitted with three neutrinos, particularly 
if one includes the Kamiokande multi-GeV results. 

The experimental indications in favor of neutrino mixing coming from the results of 
the solar and atmospheric neutrino experiments and of the LSND experiment imply that 
there are (at least) three different scales of Am^'s, about 10~^eV^, 10~^eV^ and leV^. 
We have considered all the possible schemes with four massive neutrinos which provide 
these three scales of Am^'s. 

We have shown that the results of the LSND experiment are not compatible with the 
limits obtained by all the other neutrino oscillation experiments in the case of a neutrino 
mass hierarchy (mi -C m2 -C -C m/i) and in the cases of neutrino mass spectra in 
which three masses are clustered in a group and one mass is separated from the cluster 
with a mass difference of the order of 1 eV, which corresponds to the range of Am^ 
relevant for the oscillations observed in the LSND experiment. 
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We have also shown that only two possible spectra of neutrino masses, (A) and (B) 
(see (p2D ), with two pairs of close masses separated by a mass difference of the order 
of 1 eV are compatible with the results of all neutrino oscillation experiments. If the 
neutrino mass spectrum (A) is realized in nature, the neutrino mass that is measured in 

/?-decay experiments coincide with the "LSND mass" . If the massive neutrinos are 
Majorana particles, in the case of scheme (A), the experiments on the search for (/5/3)ojy 
decay have good chances to obtain a positive result. 

Finally, we want to remark that, if the experimental indications in favor of neutrino 
oscillations are confirmed, the neutrino mass spectrum is very different from the mass 
spectra of quarks and charged leptons. This is, however, not so astonishing, because at 
least four neutrinos, one of which is sterile, are necessary in order to explain the results 
of all neutrino oscillation experiments. 

Note Added 

After this work was finished and reported at the Neutrino '96 Conference in Helsinki, 
the preprint TMUP-HEL-9605 by N. Okada and O Yasuda appeared (hepph/9606411). 
Some of our results are also contained in this paper. 
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Figure 1: Value of for the Kamiokande sub-GeV data (see Eq.(0)) as a function 
of Am?. The shadowed horizontal band represents the 90% CL limits of the double ratio 



R for the sub-GeV events determined in the Kamiokande experiment |11 . 
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Figure 2: Exclusion regions al 90% CL in the A^.g-Am^ plane for small |f/e4| and {U^^l 
in the model with mixing of four neutrinos and a mass hierarchy discussed in Section ^. 
The regions excluded by the BNL E776 and KARMEN z/^ — > i/g appearance experiments 
are bounded by the dash-dotted and dash-dot-dotted curves, respectively. The dashed 
line represents the results of the Bugey experiment. The curve passing through the circles 
is obtained from the results of the Bugey, CDHS and CCFR84 experiments using Eq.(^). 
The region allowed by the LSND experiment is shown as the shadowed region limited by 
the two solid curves. 
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